Entangled-state cycles from conditional quantum evolution 
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A system of cascaded qubits interacting via the oneway excliange of photons is studied. While for 
general operating conditions the system evolves to a superposition of Bell states (a dark state) in 
the long-time limit, under a particular resonance condition no steady state is reached within a finite 
time. We analyze the conditional quantum evolution (quantum trajectories) to characterize the 
asymptotic behavior under this resonance condition. A distinct bimodality is observed: for perfect 
qubit coupling, the system either evolves to a maximally entangled Bell state without emitting 
photons (the dark state), or executes a sustained entangled-state cycle — random switching between a 
pair of Bell states while emitting a continuous photon stream; for imperfect coupling, two entangled- 
state cycles coexist, between which a random selection is made from one quantum trajectory to 
another. 

PACS numbers: 42.50.Dv, 42.50.Lc, 42.50.Pq 



I. INTRODUCTION 

Quantum entanglement is a feature of quantum me- 
chanics that has captured much recent interest due to its 
essential role in quantum information processing 0. It 
may be characterized and manipulated independently of 
its physical realization, and it obeys a set of conserva- 
tion laws; as such, it is regarded and treated much like a 
physical resource. 

It proves useful in making quantitative predictions to 
quantify entanglement. When one has complete informa- 
tion about a bipartite system — subsystems A and B — 
the state of the system is pure and there exists a well 
established measure of entanglement — the entropy of en- 
tanglement, evaluated as the von Neumann entropy of 
the reduced density matrix. 



E{\<j>)AB)^Tr{pA\0g^PA), 



(1) 



with PA = Tr^dt/)) ((/)|^b). This measure is unity for 
the Bell states and is conserved under local operations 
and classical communication. Unfortunately, however, 
quantum systems in nature interact with their environ- 
ment; states of practical concern are therefore mixed, in 
which case the quantification of entanglement becomes 
less clear. 

Given an ensemble of pure states, with prob- 

abilities {pi}, a natural generalization of E{\(j))AB) is its 
weighted average J2iPi^i\'t'i)AB)- A difficulty arises, 
though, when one considers that a given density opera- 
tor may be decomposed in infinitely many ways, leading 
to infinitely many values for this average entanglement. 
The density operator for an equal mixture of Bell states 
= (|0)a|0)b ± \1)a\1)b)/V2, for example, is iden- 
tical to that for a mixture of |0)a|0)b and yet 
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by the above measure the two decompositions have en- 
tanglement one and zero, respectively. 

Various measures have been proposed to circumvent 
this problem, most of which evaluate a lower bound. One 
such measure, the entanglement of formation, Ep{p) 0], 
is defined as the minimal amount of entanglement re- 
quired to form the density operator p, while the entangle- 
ment of distillation, Eu{p) |3, is the guaranteed amount 
of entanglement that can be extracted from p. These 
measures satisfy the requirements for a physical entan- 
glement measure set out by Horodecki et al. 0. They 
give the value zero for pab = (|*^)($"'"| + l*I'~)(*I'~l)/2, 
which might be thought somewhat counterintuitive, since 
this state can be viewed as representing a sequence of 
random "choices" between two Bell states, both of which 
are maximally entangled. This is unavoidable, however, 
because assigning pab a non-zero value of entanglement 
would imply that entanglement can be generated by local 
operations. The problem is fundamental, steming from 
the inherent uncertainty surrounding a mixed state: the 
state provides an incomplete description of the physical 
system, and in view of the lack of knowledge a definitive 
measure of entanglement cannot be given. 

An interacting system and environment inevitably be- 
come entangled. The problem of bipartite entanglement 
for an open system is therefore one of tripartite entangle- 
ment for the system and environment. Complicating the 
situation, the state of the environment is complex and 
unknown. Conventionally, the partial trace with respect 
to the environment is taken, yielding a mixed state for 
the bipartite system. If one wishes for a more complete 
characterization of the entanglement than provided by 
the above measures, somehow the inherent uncertainty 
of the mixed state description must be removed. 

To this end, Nha and Carmichael Sj recently intro- 
duced a measure of entanglement for open systems based 
upon quantum trajectory unravelings of the open system 
dynamics 6J. Central to their approach is a considera- 
tion of the way in which information about the system 
is read, by making measurements, from the environment. 
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The evolution of the system conditioned on the measure- 
ment record is followed, and the entanglement measure is 
then contextual — dependent upon the kind of measure- 
ments made. Suppose, for example, that at some time t 
the system and environment are in the entangled state 

I0>=E^mI'^»>5|0,>£- (2) 

A partial trace with respect to E yields a mixed state 
for S. If, on the other hand, an observer makes a mea- 
surement on the environment with respect to the basis 
obtaining the "result" \(f>k)E, the reduced state 
of the system and environment is 

l-^') = l0'>s|0fc)B, (3a) 
with conditional system state 

\(t)')s =^c,^k\(f>i)s/VPk, (3b) 

i 

where pk = J2i ki.fcP is the probability of the particular 
measurement result. Thus, the system and environment 
are disentangled, so the system state is pure and its bi- 
partite entanglement is defined by the von Neumann en- 
tropy, Eq. (QJ . Nha and Carmichael jsj apply this idea to 
the continuous measurement limit, where \4>')s executes 
a conditional evolution over time. 

In this paper we follow the lead of Nha and Carmichael, 
also Carvalho et al. 7], not to compute their entangle- 
ment measure per se, but to examine the entanglement 
dynamics of a cascaded qubit system coupled through the 
oneway exchange of photons. The system considered has 
been shown to produce unconditional entangled states — 
generally a superposition of Bell states — as the steady- 
state solution to a master equation 01 . For a special 
choice of parameters (resonance) , a maximally entangled 
Bell state is achieved . . . except that the approach to the 
steady state takes place over an infinite amount of time. 

Here we analyze the conditional evolution of the qubit 
system to illuminate the dynamical creation of entan- 
glement in the general case, and to explain, in partic- 
ular, the infinitely slow approach to steady-state in the 
special case. We demonstrate that in the special case 
the conditional dynamics exhibit a distinct bimodality, 
where the approach to the Bell state is only one of two 
possibilities for the asymptotic evolution: the second we 
call an entangled-state cycle, where the qubits execute a 
sustained stochastic switching between two Bell states. 
Though involving just two qubits and elementary quan- 
tum transitions, the situation is similar to that of a bi- 
modal system in classical statistical physics in the limit 
of a vanishing transition rate between attractors. 

The physical model of the cascaded qubit system is 
presented in Sec. ^ and the quantum trajectory unrav- 
eling of its conditional dynamics in Sec. IIIII In Sec. IIVI 
we analyze the quantum trajectory equations to demon- 
strate bimodality and the existence of entangled-state cy- 
cles. Finally, a discussion and conclusions are presented 
in Sec. |V| 



II. THE CASCADED QUBIT SYSTEM 

In this section we briefly outline the physical model 
for the cascaded qubit system to be analyzed. A more 
detailed description, together with the techniques and 
assumptions used to derive the model master equation 
presented here, is available in 



A. Physical Configuration 

The system considered consists of two high-finesse op- 
tical cavities, each containing a single tightly-confined 
atom, the cavities arranged in a cascaded configuration 
with unidirectional coupling from cavity 1 to cavity 2 
(Fig. For simplicity, we consider the cavity modes 
to be identical, with resonance frequency Wcav and field 
decay rate k. Inefficiencies and losses in the coupling 
between the cavities are modeled by a real parameter e, 
< e < 1, with perfect coupling corresponding to e = 1. 
The atoms are assumed to have five relevant electronic 
levels, of which two ground states, |0) and |1), represent 
an effective two-state system, or qubit. 

For each atom, the cavity field in combination with 
auxiliary laser fields (incident from the side of the cavity) 
drives two separate resonant Raman transitions between 
states |0) and |1). An additional laser field coupled to 
the |0) ^ |t) transition provides a tunable light shift 
of the energy of state |0). All fields are assumed far 
detuned from the atomic excited states, so these states 
may be adiabatically eliminated and atomic spontaneous 
emission ignored. Under the further assumption that the 
cavity field decay rate is much larger than the transition 
rates between |0) and |1), the cavity fields may also be 



(a) 




|0> 



FIG. 1: (a) A pair of cascaded cavities. Cavity 1 and Cavity 2, 
each contain a single trapped atom; a unidirectional coupling 
between the cavities is realized by Faraday isolators F. (b) The 
atomic excitation scheme couples two stable ground states, |0) 
and to three excited states, |r), |s), and \t}. 
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adiabatically eliminated to yield a master equation for 
the reduced two-atom density matrix p, 

P^Cp = {2r,prI-r\r,p- pklk. 



1=1.2 



RiP:R\ 



R2,pR\ 



with 



(4) 



(5) 



where ai- = (|0)(l|)i, and \l3r.i\^ / k and \(3s,i\'^ /n are the 
rates of |0)i and |0)i transitions, respec- 

tively. 

By virtue of the cavity output, the system is an open 
system and solutions to master equation Q generally de- 
scribe mixed states. Under appropriate conditions, how- 
ever, the system evolves to a pure and entangled steady 
state. 



B. Steady State 

If the coupling between cavities is perfect (e = 1) and 
the parameters of the subsystems are the same (/3r,i — 
Pr,2 — Pr, /3s, 1 — Ps,2 — Ps) then the steady state is the 
pure state 



s) 



1 



(/3;|00) + /3:|11)), (6) 



where we use the abbreviated notation |00) = |0)i|0)2 
and |11) = |l)i|l)2. Then when (3^ = (3s, which we shall 
refer to as the resonance condition, the steady state is 
a maximally- entangled Bell state. This may seem to be 
ideal, but a problem arises when we consider the eigen- 
values of the operator C. Specifically, the characteristic 
time for the system to reach steady state, r — |Re(A2)|~^, 
where A2 denotes the eigenvalue of C with smallest (in 
magnitude) non-zero real part, approaches infinity as the 
resonance condition is approached. This is shown by the 
plot in Fig. 121 Thus the master equation itself, in par- 
ticular its steady state, offers limited insight into the be- 
havior of the system at resonance. We wish to learn 
more about this special case; in particular, how does the 
entanglement develop dynamically. Also, if additional 
information is factored into the description, by making 
measurements on the environment, can we better char- 
acterize the long term behavior, or possibly find perfect 
entanglement after a finite time? We demonstrate that 
quantum trajectory theory can provide answers to these 
questions. 



III. QUANTUM TRAJECTORIES 



to the environment. The first is obvious - the output 
from Cavity 2. To measure this output, let us assume 
the existence of an ideal photon detector in the path of 
the output from Cavity 2; we call it Detector 1. 

The second point of coupling to the environment is 
more subtle. Our model does not assume the inter-cavity 
coupling to be perfect; only a fraction e of the output pho- 
ton flux from Cavity 1 makes it into Cavity 2. Physically, 
this loss may be caused, for example, by non-ideal trans- 
missivity of the Faraday isolators or by absorption in the 
cavity mirrors. These imperfections cause photons to be 
scattered into the environment in some uncontrollable 
fashion. Formally, though, this is equivalent to assum- 
ing that the apparatus is ideal, except that there exists 
a beamsplitter between the cavities, as drawn schemati- 
cally in Fig. 13 We therefore further assume the existence 
of a second photon detector to collect photons reflected 
by this beamsplitter; we call it Detector 2. 

We now proceed to develop the quantum trajectory 
formalism for the cascaded qubit system. In this ap- 
proach the system is described by a pure state which is 
dependent on (conditioned on) the counting histories, or 
records, of Detectors 1 and 2. Firstly, we rewrite the 
master equation in a form suitable for translation into 
the quantum trajectory language. We reexpress Eq. Q 
in the form 



p = {Co+S)p, 



(7) 



with 



Cop = -i ^Hq, , 



-J2(c}ap + pC}C,), (8a) 

1=1,2 

(8b) 



1=1,2 




1.50 



As with any open system, the first step in unraveling 
the master equation is to identify the points of coupling 



FIG. 2: The relaxation time r = |Re(A2)| ^ plotted as a func- 
tion of |/3s//3r|. Note the singularity at resonance, |/3s//3r| = 1. 
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FIG. 3: Conceptual photon detectors, Detector 1 and Detec- 
tor 2, used for unravelling the master equation. The Faraday 
isolators are omitted for clarity. 



where 



Ho 



i?2 



V2(l-e)i?i, 
i^fe (rIRi - r\R2 



(9a) 
(9b) 
(9c) 



Then, within quantum trajectory theory, the evolution of 
the system is described by a pure state which evolves 
under the non-Hermitian effective Hamiltonian 



.1 



HcH — Hq — i- ^ CjCi, 

1=1.2 



(10) 



the continuous evolution interrupted at random times by 
quantum jumps, |0) Ci\(j)}, where the jumps occur 
with probability 



Pi{t)dt 



-dt 



(11) 



in time interval {t,t + dt). Physically, the jump opera- 
tors Ci and C2 account for the reduction of the state of 
the system, given a photon count is recorded by Detec- 
tor 1 or Detector 2, respectively. Thus, within the quan- 
tum trajectory description of the coupled cavity system, 
we consider an experiment in which ideal detectors are 
employed, such that every scattered photon is detected 
and recorded. Given the history of detector 'chcks', one 
has complete information about the system state, in the 
sense that that state is always pure; hence, although the 
solution to the master equation is generally mixed, one is 
able to characterize the entanglement in an unambiguous 
(conditional) fashion 

Consider the special case where the coupling between 
the cavities is optimal (e = 1). In this case there is only 
one output from the system, that from Cavity 2, recorded 
by Detector 1. Standard numerical algorithms 01 have 
been used to simulate typical quantum trajectories for 
various values of \(is/Pr\- Specifically, we consider the 
evolution of the conditional expectation of the operator 
product (Ti^z^2,z, where ai^z is the Pauli operator diago- 
nal in the (|0)i, |l)i)-representation. 



This expectation has a number of convenient properties; 
for example, the steady-state value 



\(^l,z<^2,. 



1, 



(13) 



regardless of the value of \PslPr\, which makes it easy 
to compare rates of convergence to the steady state for 
different system parameters. 

Figure 01 contrasts the solution to the master equa- 
tion and a single quantum trajectory. The solution to 
the master equation exhibits a completely smooth evolu- 
tion that tends asymptotically towards the steady state. 
The quantum trajectory, on the other hand, undergoes 
a sequence of switches between two extreme values of 
(<5'i.zO'2,z), which occur at each photon detection. Pro- 
vided the parameters are chosen away from resonance, 
the photon detections eventually stop and the trajectory 
settles into the steady state with {(7i^z^2,z) = 1; the 
steady state is clearly a dark state. At resonance, how- 
ever, the photon detections may continue indefinitely. 
Physically, this seems plausible, since it simply implies 
that the atoms continue to switch between states |0) and 
|1), scattering one photon with each transition. At reso- 
nance, apparently, a unique equilibrium dark state can- 
not be established. The cyclic behavior that replaces it 
is completely invisible if we consider only the ensemble 
average-a vivid demonstration of how single quantum 
trajectories can provide additional insight into the evo- 
lution of an open quantum system. 



IV. ENTANGLED-STATE CYCLES 

The oscillatory behavior featured in Fig. 0] hints at a 
simple cyclic process. In fact, it is simple enough that 
we can understand why it occurs without resorting to 
numerics. In this section we formulate a graphical de- 
scription of individual trajectories. 




ii) 1.0 



-1.0 



2 4 6 
time 



10 



ID 



ID 



.10), = -|0), 



(12) 



FIG. 4: Evolution of the ensemble average (i) compared with 
that of a single quantum trajectory (ii), for e = 1. The off- 
resonance case (Ps 7^ /3r) is shown to the left and compared 
with the resonant case (/3s = Pr) to the right. Time is mea- 
sured in units of {Pr/y/K)~ . 
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A. The Cascaded System Phase Space 

Figure2]dcmonstrates that the conditional expectation 
(<5'i,zO'2,z) is conserved during the periods of evolution 
between quantum jumps. The positively and negatively 
correlated subspaces 



{\<t>) ■■ {fyi,z»2,z) - ±1} 



(14) 



are coupled only through quantum jumps. Noting that 

=span{|00),|ll)}, span{|10), |01)} (15) 

are each 2-dimensional (assuming real amplitudes with- 
out loss of generality), we manage to break up a 4- 
dimensional space into two 2-dimensional planes, linked 
to one another by the quantum jumps. We refer to this 
representation as the cascaded system phase space. Tra- 
jectories within it can be viewed as lines moving contin- 
uosly within either plane and jumping discontinuously 
between the planes. 



B. Evolution Between Quantum Jumps 

We use phase space portraits within iJ+ and to 
characterize the behavior of the system, where for the 
sake of simplicity, and without loss of generality, we are 
assuming f3r and Ps to be real. We define 



(16) 



and scale time by setting firj \fti = 1. The master equa- 
tion then takes the form (e = 1) 



2{pR\R2 - R2pR\ + rIRip - RipRl), (17) 



where 



Ri = (7, + r&J 



(18) 



The resonance condition is now r — 1. 

It is useful to convert to a matrix notation, such that a 
pure state \(j>) of the system is represented by a 4- vector, 

10) = (cii,cio,coi,coo)^ = ^ Cy|ij), (19) 
and system operators are written as 4 x 4 matrices, e.g.. 



Ri 



and 



/ r 

r 

10 

VO 1 



i?2 



r 

10 

r 

10 



(20) 




(21) 



The evolution of 10) under H^s is written as a linear dif- 
ferential equation in four variables. 



dt' 




(22) 



As noted above, this evolution is constrained within ei- 
ther or E^ . Thus we can write |0) as a vector sum of 
two orthogonal components j^"*") G E'^ and \(j>~) G E~ , 
10) = a\(f)^) + b\(j)~), to obtain the decoupled dynamics 



-ir) = 



-2 2r 
2r -2r2 

,2\ 



(l + r2) 2r2 
2 -(l + r2) 



10- 



(23a) 
(23b) 



Eigenvectors of the two dynamical matrices correspond 
to states of the system that are preserved under the evo- 
lution between quantum jumps. Note, however, that it 
does not necessarily follow that such a state is a steady 
state of the quantum trajectory evolution as a whole; it 
must eventually experience a quantum jump if its norm 
decays — i.e., the corresponding eigenvalue is not zero. 
Recall from quantum trajectory theory that the prob- 
ability for a state not to jump prior to time t is given by 
its norm 

For the systems of equations given above we find the 
following (unnormalised) eigenstates and eigenvalues: 

(i) |0i) = |00) -hr|ll>, Ai 0; this is the steady state 
of the system for r < 1. 

(ii) 102) = -r|00)-|-|ll), A2 = -2(1-1- r2); this state in 
E'^ is orthogonal to |0i) and must eventually jump 
to a state in E^ . 

(iii) 103) = r|10) + |01), A3 = -(r - 1)2; this state in 
E^ must eventually jump to a state in _E+ unless 
r = 1; in the latter case it plays no role once an 
entangled-state cycle is initiated (see below). 

(iv) 104) ^ |10) - r|01), A4 = -(r + 1)^; this state in 
E^ must eventually jump to a state in 

In the special case of resonance, r = 1, there are two 
independent steady states, |0i) and |03), which helps to 
explain the failure of the master equation evolution to 
approach a unique steady state. It also suggests a fun- 
damental feature of the indefinite switching, the cyclic 
behavior, revealed by individual quantum trajectories: 
during such an entangled- state cycle, the system state 
must remain orthogonal to |0i) and |03). We verify this 
shortly, after examining the trajectory evolution away 
from resonance, where the steady state |0i) is always 
reached for perfect inter-cavity coupling. 
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C. Quantum Trajectories for r < 1 

Typical quantum trajectories for r = 0.5 are shown in 
Figs. Eland El where the and subspaces are drawn 
as circular planes. Normalized states are located on the 
circumferences of the circles. The Bell states 



= (|00)±|11))/V2, 
= (|01)±|10))/V2, 



(24a) 
(24b) 



lie at intersections of the circumference with the dotted 
lines as shown. 

Between quantum jumps, under the influence of the 
non-Hermitian Hamiltonian Hos, the norm of the state 
decays and the point representing it within the phase 
space moves to the interior of one of the circles. Quantum 
jumps cause a switch from to E~ or vice-versa. They 
are represented by the lines connecting the two planes, 
where for illustrative purposes, the system state is renor- 
malized after each quantum jump; thus jumps terminate 
at points on the circumference of the circles. 

We restrict ourselves to separable initial states located 
in one or other of the two subspaces; for example, the 
states 1 00) and |10), respectively, are considered in Figs. El 
and El 



1. Effect of quantum jumps 

The action of the jump operator Ci on states located 
in E'^ (with renormalization) is 





signj cii - rcpo} / -1 

V2 1 




while the action of Ci on states in E is 




(25a) 




signjcio - cqi} 






V 1 



(25b) 



Thus, when a quantum jump occurs, any state within 
£■+ collapses onto the Bell state \4>) = ij^f") in E^ , 
while any state within E~ collapses onto the state |0) — 
±(|00) - r|ll))/Vl + ^2 in E+. 



while with probability = 1 - P+^- = 

the system evolves to the steady state without any 
photon emissions. 

If a jump from to \(^) = ±(|00) - r|ll))/\/l + r-^ 
has just occurred, then by the same argument one shows 
that the probability of a future quantum jump to E^ is 
4r2/ (l-|-r'^)2, or, alternatively, the probability of reaching 
the steady state after such a jump is 1 — 4r2/ (1 + r^)^ = 
[(l-r2)/(l + r2)]2. 

Consider now an initial state in E~ ^ |</)~(0)) = cj^a) + 
d|04), for some (real) {c, c?}. Owing to the instability of 
both 103) and |04) for r < 1, an eventual quantum jump 
is guaranteed; thus. 



1. 



(27) 



Armed with this information, we move to an explanation 
of the quantum trajectories displayed in Figs. El and El 



3. Initial states jOO) and jlO) 

In Fig. El we plot three typical phase-space trajectories 
for r = 0.5 and \(j)+{0)) = |00). Fig. Ela) illustrates 
the case where the system evolves directly to the steady 
state \4>i)- The probability of this event is |(OO|0i)|^ = 
l/{l + r^) = 0.8, so it is the most likely occurrence for the 
chosen parameters. If a first quantum jump does occur, 
then typical trajectories are shown in Figs.E^b) and (c). 
Following the jump to \(j)) = — j^*") in E~ , a second jump 
returning the state to i?+ is guaranteed. For r = 0.5, this 
leaves the system in the state 10) = 0.89|00) - 0.45|11), 
from which the probability of a further cycle of jumps is 
4r2/(l -I- r^)^ = 0.64. Thus, after a first quantum jump 
cycle, it is most likely that further cycles will follow, as 
seen in Figs. EJb) and (c), where in both cases a total 
of five cycles (ten photon detections) occur before the 
system finally reaches the steady state. 

In Fig. El we plot three typical phase-space trajectories 
for r — 0.5 and \4>~{0)) — |10). In this case, at least 
one quantum jump is certain to occur, following which 
the probability of further jumps is Ar"^/ (1 -l-r^)^ = 0.64, 
as above. So for this initial condition, the most likely 
outcome is a sequence of quantum jump cycles following 
a first guaranteed photon detection. In Fig. El (a) only 
the first detection occurs, while in Figs.E^b) and (c) this 
detection is followed by a sequence of cycles before the 
steady state is eventually achieved. 



2. Jump probabilities 

Consider an initial normalized state in E^ , |0+(O)) = 
a\(j)i) + b\(j)2), for some (real) {a,b}. Given that |0i) is 
a steady state of the evolution between quantum jumps, 
the probability of an eventual quantum jump to E~ is 

^+^- = |(0+(O)|02)|' = 6^ (26) 



D. Quantum Trajectories for r = 1 

The case r = 1 is of particular interest. The normalized 
eigenstates of the evolution between quantum jumps are 
the Bell states |0i) = |$+), l^a) = |$-), \h) = l*+), 
and 104) — \'^~)- The eigenvalues are Ai = A3 = and 
A2 = A4 = —4. The action of the jump operator Ci on 
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<^I11>,<^I10> <^lll>,<^IIO> 




<^I11>,<^I10> <^I11>,<$I10> 




<(^lll>, <^IIO> <(^lll>, <^IIO> 



FIG. 5: Examples of phase space trajectories for r = 0.5 with FIG. 6: Examples of phase space trajectories for r = 0.5 with 
initial state |</)^(0)) = |00) and e = 1. See text for description. initial state \4>~ {Q)) = |10) and e = 1. See text for description. 



states within simplifies to 

Ci M -.sign{cn-coo}|*-), (28a) 
V Coo / 

and its action on states within E~ to 

Ci y ^sign{cio-coi}|$-). (28b) 



For r = 1, photon detections, if they occur, are associated 
with collapses onto one of two maximally-entangled Bell 
states. 

For initial states |(/)+(0)) and |0"(O)) in E+ and E' , 
respectively, the system evolves continuously, without the 
emission of any photons, to — l^"*") and l^a) — \'^~^), 
with probabilities |(0+(O)|$+)|^ and |((/)-(0)|^'+)|^ Al- 
ternatively, a photon is detected with associated quan- 
tum jump to l^*") in E~ or |$~) in E~^. In this case, 
as both terminal states are unstable under the between- 
jump evolution, a second detection and quantum jump 
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must follow. According to Eqs. (|28a|l and Ij28b|l this sim- 
ply exchanges l^*^) for 1$^) and vice-versa. Hence, a 
perpetual switching between Bell states l^*^) and |(f>^) 
occurs. We designate this behavior an entangled- state 
cycle. 

Thus, at resonance we find a distinctly bimodal be- 
havior. The system either evolves into a maximally- 
entangled Bell state without emitting photons, or an 
entangled-state cycle is initiated under which the system 
switches indefinitely between orthogonal Bell states while 
emitting a continual stream of photons. As an aside, 
such behavior can be regarded as a quantum measure- 
ment that distinguishes the Bell states from l^*^). 

The two alternative outcomes of the quantum trajec- 
tory evolution are illustrated in Figs. [7\ and |S1 for the 
initial states |(/)+(0)) = |00) in E+ and |0"(O)) = |10) in 
E~ , respectively. With this choice of initial states there 
are equal probabilities for reaching the steady states, 
1*^) [Fig-CJa)] and 1*+) [Fig.|HIa)], and for commencing 
an entangled-state cycle [Figs.Cfb) andlH^b)]. Note that 
once an entangled-state cycle is initiated, the trajectory 
remains in a plane orthogonal to the lines defining |$^) 
and 1^'"*'); the cycle continues indefinitely 
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FIG. 7: Examples of phase space trajectories for r = 1 with 
initial state |(j!)^(0)) = |00) and e = 1. See text for description. 



E. Imperfect Intercavity Coupling 

Our original model allowed for the possibility of imper- 
fect intercavity coupling, through the parameter e and 
the jump operator C2 which describe the effects of pho- 
ton loss in propagation between the two cavities. Fo- 
cusing on the resonant case (r — 1), we now consider 
the situation in which e < 1. Typical trajectories for 
e = 0.5 are shown in Figs. |5fa) and lior a'l. with the two 
photon count records shown in frames (b) and (c) of the 
figures. Remarkably, entangled-state cycles persist, but 
now the system settles into one or other of two distinct 
cycles, involving either the symmetric or antisymmetric 
Bell states. 

To understand the behavior, consider the forms of the 
operators involved; in particular, for r = 1, we have ef- 
fective Hamiltonian 
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FIG. 8: Examples of phase space trajectories for r = 1 with 
initial state |(;/>~(0)) — |10) and e = 1. See text for description. 
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FIG. 9: (a) Example of a phase space trajectory for r = 1 
with initial state 10^(0)) — |00) and e = 0.5; the system 
eventually settles into the |<1?^) <^ l^* ) entangled-state cycle, 
(b) and (c) Photon counts for Detectors 1 and 2, respectively. 
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FIG. 10: (a) Example of a phase space trajectory for r = 1 
with initial state \(f>'^{0)) — |00) and e — 0.5; the system 
eventually settles into the \$^) ^ I'l'^) entangled-state cycle, 
(b) and (c) Photon counts for Detector 1 and 2, respectively. 
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Significantly, these operators commute with one another, 

[Ci, (72] = [Cl, -ffeff] = [C2, Hes] = 0. (31) 

Their operation upon the Bell states is given by 

-iHe«\<^^) = -2(1tV^)|*^) = A±|$±), (32a) 
-z^eff|*±) = -2(lT\/^) 1*^) =A±|*±), (32b) 



Cil**) 



= ±v^(i^|*±) 
= ±V2(T^|$±) 



(33a) 
(33b) 
(33c) 
(33d) 



Thus, the Bell states are eigenstates of HeS, and the jump 
operators interchange Bell states in and E~: each 
jump operator converts the symmetric (antisymmetric) 
Bell state in to the symmetric (antisymmetric) Bell 
state in E~ and vice-versa. 

Now, let us consider a particular quantum trajectory 
for which a total of n jumps occur, separated by the time 
intervals {Ati : i — 1, ... ,71}. For an initial state |(/)o), 
the (unnormalized) state at the conclusion of the n jumps 
is written as 



t/716 



iH^„At„ 



Jie 



iH^ffAt 



where each Jj is either Ci or 6*2. Since all operators in 
the string acting on lipg) commute, this expression can 
be rewritten in a variety of forms, two of which prove to 
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be especially useful in explaining the distinct behaviors 
illustrated by Figs. 1^ and 1101 In the first case, we may 
write 

mii)^CTe-^"^"'{ci\M), (34a) 

passing all I ocurrences of Ci to the right and all m oc- 
currences of C2 to the left {I + m = n); in the second we 
write 

\^t)in)=CiCT{e-^"^-'\^o)), (34b) 

where all jump operators are passed to the left. 

The arbitrary (pure) initial state can be expressed as 
a superposition of Bell states, 

|0o) = a|*+) + b\^-) + c\^+) + (35) 

where a, b, c, and d are expansion coefficients, generally 
complex. Substituting this expansion into Eqs. (|34a|l and 
(|34b|l ■ and using Eqs. (|32a|l - lj33d|) — assuming for simplic- 
ity that I and m are even — the two forms for the state 
|0t) are 

(i) cx e-'^-*|0o)', (36a) 
« C'M, (36b) 

where 

\M' = AV(a|$+)+c|vI/+)) + AL(6|<&-)+d|*-)), 

(37a) 

\(t>ty = e^+*(a|$+)+c|*+)) + e^-*(&|$-)+d|*-)). 

(37b) 

Observe now that the ratio of the eigenvalues satisfies 

A+/A_ = (1-V^)/(1 + V^)<1. (38) 

It follows that \(f>t){i) and \(f>t)(n) allow us to predict quite 
distinct asymptotic behaviors for the system state. For 
sufficiently large I, the contribution to \4't){i) from the 
symmetric Bell states is negligible compared with the 
contribution from the antisymmetric Bell states, in which 
case, using Eqs. (|36a|l and (|37a|l . 

|</'t)(i)~e^-*Ai.(6|$-)+d|*-)). (39) 

The system is locked into a cycle between the two anti- 
symmetric Bell states, the situation illustrated in Fig. |5| 
(for e — 0.5, |A+/A_| = 0.17). In contrast, for sufficiently 
large t, the contribution to \4>t)(u) from the antisymmet- 
ric Bell states is negligible compared with that from the 
symmetric Bell states, and using Eqs. (j36bp and (|37b|) . 

|0,)(ii)-e^+*(^^y(fe|a>+>+d|*+». (40) 

The system is locked into a cycle between the two sym- 
metric Bell states, as shown in Fig. ^| 



Which of the two cycles is chosen in a particular re- 
alization of the photon counting record is random, as 
is the time taken to settle into the cycle. Effectively, 
the decision is the outcome of a competition between 
the periods of evolution between quantum jumps and 
the jumps themselves — specifically, those associated with 
photon counts at Detector 1. Considering Eqs. H37a|) 
and (|38|) . we see that every count at Detector 1 results 
in an increased probability to find the system in one of 
the antisymmetric Bell states. On the other hand, from 
Eqs. I|37b|l and H38|l . the periods of evolution between 
counts have the reverse effect — they increase the prob- 
ability for the system to be found in a symmetric Bell 
state. The critical factor that decides which tendency 
wins is the number of photon counts occuring at Detec- 
tor 1 over a given (substantial) interval of time. If there 
are many, as in Fig. Efb) , the entangled-state cycle be- 
tween antisymmetric Bell states wins out; if there are few. 
Fig. Iior b). the cycle between symmetric Bell states oc- 
curs. The same decision mechanism is observed in other 
examples Note that counts at Detector 2 are not 

involved — not directly at least. They do figure indirectly 
as a mechanism reducing the average number of counts 
at Detector 1; indeed, they are the ultimate source of the 
asymmetry reflected in the ratio A+/A_ < 1. 

As the system approaches a particular cycle the quan- 
tum trajectory evolution tends to reinforce the estab- 
lishment of the cycle. Close to the antisymmetric cycle, 
the evolution between jumps is dominantly governed by 
A_ = —2(1-1- ^/e) and is therefore relatively fast. This 
leads to frequent photon counts at Detector 1 [Fig. Elb)] . 
Close to the symmetric cycle, the between-jump evolu- 
tion is dominantly governed by A+ = —2(1 — y/e), hence 
is relatively slow. Photon counts at Detector 1 become 
much less frequent [Fig. lTUT b)]. 

From the dramatic difference in count rates at Detec- 
tor 1 for the two cycles, it is clear that one can deter- 
mine which entanglement cycle the system evolves to for 
a particular realization. However, without knowledge of 
the record of photon counts at Detector 2, which by def- 
inition we do not have, one can not know where on the 
cycle the system is, i.e., whether the state is in or 
. Thus, the ensemble average state of the system is 
mixed, described by one of the density operators 

p± = ^(|$±)($±| + |*±)(vl,±|). (41) 



V. DISCUSSION AND CONCLUSIONS 

Consider a thought experiment where the cascaded 
qubit system, set to resonance, evolves freely and its en- 
tire output is collected and stored inside a black box. At 
some time the lasers driving the Raman transitions are 
turned off, so the evolution ceases. The box and qubits 
are separated and moved to causally disconnected regions 
of space time. Let Alice and Bob be standard observers 
of the qubits, and give Eve jurisdiction over the box. 
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We can now ask, how much entanglement exists be- 
tween the qubits of AUce and Bob? While this is simply 
a roundabout way of asking how entanglement evolves, 
it helps elucidate some of the key concepts behind the 
quantum trajectory measure of entanglement. Conven- 
tional entanglement measures are based upon an analysis 
of the density matrix at this time. They throw away the 
box and look at the system of qubits alone — they disre- 
gard Eve and view the system from the perspective of 
Alice and Bob. 

Yet in general every interaction between two objects 
entangles them, and as the qubit system and box inter- 
acted in the past, their states are intertwined. Neither 
possess an independent reality, and neither, considered 
alone, can be completely described. Eve's box contains 
information, which, if discarded, adds entropy to the 
qubit system of Alice and Bob. This entropy is the source 
of ambiguity in the quantification of entanglement. From 
this point of view, as noted in the introduction, the prob- 
lem of bi-partite entanglement in an open system relates 
to that of tri-partite entanglement in a closed one. To 
completely characterize the entanglement of the present 
example, in addition to the entanglement between Alice 
and Bob, we must consider their entanglement with Eve. 

A quantum description of the box is impractical, but 
it is feasible to extract classical information about what 
it contains, through measurement. Quantum trajecto- 
ries facilitate this, and allow us not to discard the box 
completely. In turn, the system state retains its purity, 
conditional on the classical information extracted from 
the box. With this extra information, we can extract 
more entanglement from the cascaded qubit system. 

Working from the master equation for the cascaded 
system previously it was assumed that the system 
evolved gradually into a pure state, whereby entangle- 
ment was generated. The behaviour at resonance, how- 
ever, was unclear, since there the master equation had 
two zero eigenvalues and no well-defined steady state. 
By considering the conditional evolution we have shown 
that, at resonance, asymptotically the system is either in 
the Bell state or oscillating (stochastically switch- 
ing) between two Bell states, |$~) and |^^). 

From the density matrix point of view, the latter is 
an equal mixture of Bell States and would yield no en- 
tanglement under any mixed state measure; physically, 



Alice and Bob, without collaboration from Eve, cannot 
extract any entanglement from their qubits. Suppose, 
however, that Eve opens her box to count the number of 
photons inside. Seeing whether the count is even or odd, 
she is able to deduce exactly which Bell state Alice and 
Bob's system is in. Thus, her measurement unravels the 
density operator, creating entanglement, despite the fact 
that the measurement is not causally connected to Alice 
and Bob's qubits. 

It is tempting to say that the entanglement was al- 
ways there, as a matter of fact, until one realizes that 
there are many other ways in which Eve could choose 
to measure her state, each producing a different unrav- 
elling of the qubit system and yielding a different value 
of entanglement. The entanglement facilitated by Eve's 
measurements is contextual in this sense. 

This thought experiment demonstrates why any at- 
tempt to quantify the entanglement of an open system 
from the density operator alone cannot be considered 
complete. The density operator should not be treated 
as a fundamental object, as it does not provide a com- 
plete description of the physical state. We have presented 
a simple example where oscillations between maximally 
entangled states are hidden within a separable density 
operator. The fact that the density operator contains 
entropy, implies that information about its entanglement 
with an external system was discarded at some time. In 
studying such a mixed state, there is benefit from con- 
sidering, not only the mixed state itself, but the process 
through which it was generated, and the access this po- 
tentially gives to a conditional dynamics. 

The results of this paper could be extended by employ- 
ing quantum trajectories in a broader sense. In cases 
where the results of environmental interactions cannot 
be measured, such as coupling loss, Wiseman and Vac- 
caro |lOl | have shown that only certain unravelings can 
be physically realized. A conceivable measure of entan- 
glement would take the minimum of all physically re- 
alizable unravelings. Alternatively, one might take the 
maximum of all physically realizable unravellings, which 
would measure the maximum distillable entanglement 
when local measurements on the environment are taken 
into account. 

This work was supported by the Marsden Fund of the 
RSNZ. 
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